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Abstract. We present results on the total mass and tem- 
perature determination using two samples of clusters of 
galaxies. One sample is constructed with emphasis on the 
completeness of the sample, while the advantage of the 
other is the use of the temperature profiles, derived with 
ASCA. We obtain remarkably similar fits to the M — T 
relation for both samples, with the normalization and the 
slope significantly different from both prediction of self- 
similar collapse and hydrodynamical simulations. We dis- 
cuss the origin of these discrepancies and also combine the 
X-ray mass with velocity dispersion measurements to pro- 
vide a comparison with high-resolution dark matter sim- 
ulations. Finally, we discuss the importance of a cluster 
formation epoch in the observed M — T relation. 
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1. Introduction 

There is increasing interest in using galaxy clusters char- 
acterized by X-ray observations as probes for cosmic struc- 
ture and its evolution. In one of the applications, Press- 
Schechter models (Press & Schechter 1974) are used to 
constrain the shape and amplitude of the primordial den- 
sity fluctuation spectrum at present day scales of 5 — lOh^ 1 
Mpc (e.g. Henry & Arnaud 1991, Oukbir & Blanchard 
1992, Henry et al. 1992, Eke et al. 1996, 1998, Borgani 
et al. 1999). Another important measurement is that of the 
spatial correlations of X-ray clusters or, alternatively, the 
measurement of the power spectrum of the density fluc- 
tuations in the cluster density distribution (e.g. Nichol 
et al. 1992, Romer et al. 1994, Retzlaff et al. 1998, 
Miller & Batuski 2000, Schuecker et al. 2000, Collins 
et al. 2000). In both cases, the masses of the clusters have 
to be known for a large test sample in order to apply it 
to the structure analysis. The mass determination is still 
only feasible for a limited number of well-observed galaxy 
clusters, however. Therefore, in the above approaches the 
empirically-found correlation of the mass with the more 
easily observable parameters, X-ray luminosity or intra- 



cluster gas temperature, are used to make the connection 
between theoretical modeling and observations. The mass- 
temperature relation plays a crucial role in this analysis, 
since it was found to be particularly tight. Since the time 
of the first X-ray temperature measurements it was ob- 
served that the X-ray temperature is closely connected to 
the velocity dispersion of galaxies in clusters, indicating 
that both cluster components trace the cluster potential 
in a similar way (e.g. Mushotzky et al. 1978, Mushotzky 
1984). 

The validity of this approach was investigated in nu- 
merical calculations based on A-body dynamic and hy- 
drodynamic simulations making predictions on the mass- 
temperature relation and its dispersion, as measured in 
X-ray studies (Evrard, Metzler, Navarro 1996, Evrard 

1997) . A close correlation with a dependence of M cx T 3 / 2 
was predicted. It was found in this work and also by 
e.g. Schindler (1996) that the X-ray mass determina- 
tion should be reliable (with an uncertainty in the range 
14 - 29%). Evrard et al. (1996 and Evrard 1997) argue, 
however, that the predicted mass-temperature relation has 
such a small dispersion that a mass estimate based on 
the temperature measurement only is more accurate than 
that determined on the basis of the additional knowledge 
of the gas density profile as obtained from /3-model sur- 
face brightness profile fits. These predictions from simula- 
tions have been tested recently by means of observational 
results by several authors (e.g. Bohringer 1995, Hjorth 
et al. 1998, Horner et al. 1999, Nevalainen et al. 2000). 
Disagreements have been found concerning the predicted 
slope (see e.g. Horner et al. 1999, Nevalainen et al. 2000) 
and the normalization of the mass-temperature relation. 

At the time of writing, this discrepancy has not been 
solved, and various aspects of the data analysis and 
simulations have been discussed: the correct measure- 
ment of temperature gradients in the intracluster gas 
(e.g. Markevitch et al. 1998; Irwin, Bregman, Evrard 
1999, White 2000), influence of cooling flows (Allen 

1998) , non-thermal pressure support (e.g. Schindler 1996; 
Fukazawa et al. 2000), influence of the heating of the 
intracluster medium by supernovae (e.g. Balogh, Babul, 
Patton 1999; Loewenstein 2000; Finoguenov, Arnaud & 
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David 2000, hereafter FAD), density and velocity bias in 
the galactic component of the cluster (Colin et al. 1999) 
and effects of numerical resolution in the simulations 
(Nevalainen et al. 2000). 

In this paper we reinvestigate the observational mass- 
temperature relation using an extended sample of clus- 
ters with measured temperature profiles covering a wide 
range of systems with luminosity-averaged temperatures 
from below 1 up to 10 keV. In particular, we include 22 
low temperature systems, thus improving on sampling the 
low- mass part of the M — T relation, compared to previous 
studies (e.g. Horner et al. 1999; Nevalainen et al. 2000). 
To study in detail the parameters that influence the nor- 
malization and the slope of the M — T relation, we rederive 
the M — T relation assuming an isothermal temperature 
distribution and also compare with the M — T relation for 
HIFLUGCS, a statistically complete sample of clusters, 
described in more detail below, where the isothermality 
assumption was the only choice. This sample is included 
in our study primarily to demonstrate that the result of 
the M — T relation is not affected by any selection bias. In 
addition, we discuss implications from the results of the 
M — T relation study on the epoch (redshift) of cluster 
formation. 

The paper is organized as follows. In section ^ we de- 
scribe the sample compilation, in ^ we compare the results 
obtained for the two different samples, in ^ we combine 
the X-ray mass with velocity dispersion measurements 
to provide a comparison with high-resolution dark mat- 
ter simulations and in || we discuss a correction for the 
redshift of cluster formation. Unless noted otherwise, we 
assume 57m = 1, A = and H Q = 50 km s -1 Mpc -1 , 
Pcrit,o = 4.6975 x 10 -30 g cm~ 3 throughout the paper. 

2. Data Compilation and Results 

In this paper we compare the M — T relation of two clus- 
ter samples. The first is a statistically complete sample 
of the X-ray brightest clusters, the Highest X-ray FLUx 
Galaxy Cluster Sample (HIFLUGCS), for which the selec- 
tion criteria and thus possible bias effects are well defined. 
For this sample, masses have been determined on the as- 
sumption of isothermality of the X-ray emitting plasma. 
The second sample comprises those clusters for which we 
have temperature profiles from ASCA observations and 
for which a refined mass determination is possible. For 
the analysis of both samples, we use the same definition 
of the "total mass" as the mass enclosed by the radius, 
rsoo, inside which the mean cluster mass density is 500 
times higher than the critical density of the Universe. Ac- 
counting for the observed redshift means that the critical 
density, used to calculate the overdensity, is scaled for ev- 
ery cluster according to its redshift. Later, in fitting the 
M — T relation, we correct the temperature of the clus- 
ter by dividing by (1 + z). These corrections mean that 
clusters form at different epochs at similar overdensity, 



so M3/(4irR 3 ) = 180p crit , o (l + zf. Therefore, for a fixed 
mass R oc (1 + z)' 1 , and T oc M/R oc (1 + z). This 
is strictly correct only for Om = 1. However, correction 
for the observed redshift assuming lower values of £Im has 
even smaller effect on the derived parameters of the M — T 
relation, so our correction, assuming 17m = 1, can be con- 
sidered as an extreme case. As we will show below, even 
in this case, there is no difference in the results obtained 
for both cluster samples considered here and thus no dif- 
ference in the derived parameters of the M — T relation is 
expected for any value of £Im- 

2.1. HIFLUGCS 

Candidates for HIFLUGCS have been selected from re- 
cent cluster catalogs based on the ROSAT All-Sky Survey. 
They have been reanalyzed homogeneously to construct a 
complete X-ray flux-limited sample of the brightest clus- 
ters in the sky, comprising 63 clusters. Details of the sam- 
ple construction are described in Reiprich et al. (2001). 
The cluster masses have been calculated by assuming hy- 
drostatic equilibrium and isothermality of the intraclus- 
ter gas and determining the gas density profile with the 
(3 model (Cavaliere and Fusco-Femiano 1976, Gorenstein 
et al. 1978, Jones and Forman 1984). Overall cluster tem- 
peratures have been compiled from the literature, giv- 
ing preference to temperatures measured with the ASCA 
satellite. The cluster masses, M500 and the mass errors, 
determined by adding the temperature error and the er- 
rors of the fit parameter values ((3 and core radius), are 
tabulated in Reiprich et al. (2001). Since HIFLUGCS is 
purely flux-limited, we avoid any bias that could be intro- 
duced in a subjectively compiled sample, e.g. based on 
the most preferred targets selected for deep observations. 
To maintain the completeness of the sample, out of 63 
systems in HIFLUGCS, for two we have to use estimated 
temperatures based on the L — T relation by Markevitch 
(1998). We have checked, however, that excluding these 
two systems from the fit does not change the derived pa- 
rameters of the M — T relation. 

In addition, we will use also a larger sample (en- 
larged HIFLUGCS, 88 clusters), which is not strictly flux- 
limited, but contains only clusters with measured temper- 
atures. 

A linear regression is performed in log(T)-log(AT) 
space. The method allows for errors in both variables 
and intrinsic scatter (Akritas & Bershady 1996). We use 
the bisector method to determine the best-fit parameters 
(Isobe et al. 1990). Errors are transformed into log space 
by Alog(x) = log(e) (x + — x~)/(2x) where x + and x~ 
denote the upper and lower boundary of the quantity's 
error range, respectively. The values of the fit parameters 
are given m Table H 

The results, obtained for the flux-limited and the en- 
larged sample agree within the uncertainty of the fit. A 
correction for the observed redshift does not result in any 
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Table 1. Fit parameter values for the M — T relation 7 . 



sample slope 


norm, 1 keV 
10 13 M 


HIFLUGCS 


flux-limited 1.676 ± 0.054 


o ro+0.33 
°- do -0.30 


enlarged 1.636 ± 0.044 


a 74 +0.29 
J - '^-0.27 


flux-limited, p crit = p<ait(z) 1-679 ± 0.054 


q ro+0.33 
°- oz -0.31 


enlarged, p cr i t = /Ocrit(z) 1.636 ± 0.046 


a 74 +0.29 
°- '^-0.27 


Sample with fcT-profiles 


entire sample 1.781q'q9 


2 fi1+ - 38 

z - ul -0.33 


entire sample, p crit = p C rit(z) l-78io!o9 


2 64+ 39 


M > 5 x 10 13 M Q 1.58±g;gf 


q r 7 +0.41 
J ' d ' -0.35 


p > o.4 i-58j:S:Sj 


q ry+0.27 
J -°' -0.26 


Implying isothermality to sample with fcT-profiles 


entire sample 1.89j]Q" Qg 


2 4n+ - 37 

z -^ d -0.32 


M > 5 x 10 13 M Q l-74±S;ge 


3 04 +o.29 

°- uq; -0.29 


/8 > 0.4 l-66±°;°t 


q Kn+ - 21 
O.OU_ 23 



7 Errors are given at the 68% confidence level. A bootstrap 
method is used to estimate the errors. 



change in the derived parameters. The derived slope of 
the M — T relation is steeper than the value of 1.5, ex- 
pected from self-similar scaling relations. The normaliza- 
tion obtained in the hydrodynamical/A^-body simulations 
of Evrard et al. (1996) is higher than found here. No break 
in the M — T relation is visible over the whole range of 
temperatures. 

Note, however, that the determination of the mass 
value itself depends on the temperature, therefore some 
care has to be taken in the interpretation of the fit re- 
sults. In the discussion, which follows below, we take this 
effect into account. 



2.2. The sample with known temperature gradients. 

To derive the total masses for the clusters in this sam- 
ple, we use the spatially resolved temperature profiles 
found in ASCA measurements (Markevitch et al. 1998, 
Finoguenov & Ponman 1999, Finoguenov, David & Pon- 
man 2000, hereafter FDP, FAD, Finoguenov, Jones & 
Bohringer 2001). The sample totals 39 systems with tem- 
peratures from 0.7 keV to 10 keV. This corresponds to 
a factor of 100 difference in total mass, determined at a 
given overdensity. A major difference of this sample, com- 
pared to studies of Horner et al. (1999) and Nevalainen 
et al. (2000), is an inclusion of 22 systems with tem- 
peratures spanning the range from 0.7 keV to 3.5 keV. 
This sample is therefore well suited to study the possi- 
ble break in the M — T relation, suggested by Nevalainen 
et al. (2000). 



In calculating the total cluster masses we used poly- 
tropic indices to describe the temperature profiles, omit- 
ting the cluster core, where effects of cooling may be im- 
portant. For the total mass estimates, we used the fits to 
the surface brightness profiles from ROSAT PSPC data on 
the outskirts of the clusters from Vikhlinin et al. (1999), 
FDP, FAD, Finoguenov, Jones, Bohringer (2001), thus 
avoiding the cooling zone of the cluster. We estimate the 
uncertainty of surface brightness profile fitting on the mass 
estimation as 4% and propagate this error to the total 
mass. The uncertainties in the total mass estimations are 
much larger and are due to the uncertainty in temperature 
estimates and temperature gradients. 

There could be a systematic effect in the analysis of 
ROSAT surface brightness profiles of groups, caused by 
variation of the ROSAT countrate-to-emission measure 
conversion factor at low emission temperatures (Bohringer 
1996; Finoguenov et al. 1999). To evaluate this effect, 
we compare our fits to ROSAT surface brightness profiles 
with the spectral normalizations, derived in our 3d model- 
ing of ASCA data. The latter accounts for both tempera- 
ture and metallicity effects, but has larger systematics due 
to the complex PSF. From a good agreement, found in this 
comparison, a possible systematics in the total mass cal- 
culation, resulting from our approach to measuring the (3 
values, is constrained to be less than 10%. 

Under the assumption of hydrostatic equilibrium, 
when the density profile is described by a /3-model and the 
temperature distribution is expressed in polytropic form 
(T(r) oc n^as), the total mass within the radius r = xr c 
is simply 



M tot {< r) = 3.70 x 10 ld M Q T(r)r 



3f3jx 2 



1 



(1) 



We employed Eq.[l] for our mass measurements. In 
Table | we list our mass determinations at overdensity 
500. Columns denote system (1), redshift (2), emission 
weighted temperature (3), temperature at r§oo (4), total 
mass within 7- 50 o in 10 14 M Q (5), measured r 50 o in Mpc 
(6), beta and core radius (7-8). Temperature gradients are 
tabulated in col. (9), expressed as polytropic indices (7). 
In col. (10) we cite the outer radii included in the analysis 
of ASCA data. All errors in this table are given at the 
68% confidence level. 

In calculating the emission-averaged temperatures we 
removed the effects of cooling flows and emission lines. 
Thus, these temperatures are not subject to effects dis- 
cussed in Mathiesen & Evrard (2000, hereafter ME00). 
The deviation of the measured M — T relation relative to 
the simulated one could be characterized by observing a 
higher temperature for a given mass. If a distinction exists 
between the spectral temperature measurements and the 
mass-averaged temperatures, as discussed in ME00, the 
above-mentioned discrepancy should only increase. This, 
however, is not true in the case of decreasing temperature 
profiles. We will return to this issue below. 
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Fig. 1. M — T relation for the enlarged HIFLUGCS (filled circles indicate the data with solid line indicating the best-fit). 
For comparison, the fit to the M — T relation for the flux-limited HIFLUGCS is also shown (dot-dashed line). The dashed 
line shows the result of simulations by Evrard et al. (f996). 
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Fig. 2. M — T relation (analog to Fig|l|) for the sample with temperature profiles. Crosses represent the mass determinations 
using ASCA temperature profiles and ROSAT surface brightness profile fitting. The dotted line denotes the best fit using the 
total sample, while the solid line denotes the best fit, when groups (M500 < 5 x 10 13 Mg) are excluded from the fitting. The 
dashed line shows the result of simulations by Evrard et al. (f996). The long-dashed line shows the fit for the low-temperature 
end of the sample (kT ew < 4.5 keV). The grey line shows the effect of SN preheating on the M — T relation, discussed in FAD. 
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2.3. M — T relation from spatially resolved temperatures. 

In this section we investigate how the resulting parameters 
describing the M — T relation depend on the selection of 
the sample, with the most important results listed in Table 

0- 

Without the correction for the observed redshift the 
fit to the M — T relation using the bisector method by 
Akritas & Bershady (1996) gives 

M 500 = (2.6lt° : g)10 13 x fcTeJ- 09 . 

All the errors for the M — T relation reported in this 
section were calculated from 10000 bootstrap realizations. 
The mass is expressed in Mq and the temperature in keV. 



0.14, 90% confidence interval is cited) at an overdensity of 
1000. 

Horner ct al. (1999) find a flatter M — T relation, con- 
sistent with a value of 1.5, but their sample lacks groups. 
Nevalainen et al. (2000) suggest on the basis of their 
comparison of the M — T relation, derived for hot clus- 
ters and adding a few groups, that there might be a break 
in the M — T relation, occurring below 4 keV. Since our 
data uniquely sample a temperature range from groups 
to clusters of galaxies, we are in a position to check this 
suggestion. 

The M — T relation derived without groups (M 50 o > 
5 x 10 13 Mq) is 

M 500 = (3.57l°ji)10 13 x kTt ,- om . 
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M(kT=lkeV) 10 13 M 
Fig. 3. Parameters of the bisector fit to the M—T relation. 
Contours, drawn at 1,2 and 3 a confidence level, denote 
the parameters derived from the sample with tempera- 
ture profiles. In grey we show the fit to all the data and in 
black those excluding the groups (M^qq < 5 x 10 13 Mq). 
The black point shows the value obtained in simulations 
of Evrard et al. (1996), while the black crosses denote 
the values obtained for the enlarged and flux-limited HI- 
FLU GC5 sample with the error bars shown at the 68% 
confidence intervals. 



Correcting for the redshift we obtain 



1.78 



.-.III] 



(2.64^)10 13 x kT ew 



+ 0.10 
0.09 



Thus, there is no significant effect due to the redshift 
correction for our sample, composed from nearby objects, 
and we will omit it in further relations. 

Our results on the slope of the M — T relation are in 
agreement with findings of Nevalainen et al. (2000), who 
determined a significantly steeper slope than 1.5 (1.79 ± 





T relation, 



M(kT=3.5 keV) 10 14 M 

Fig. 4. Parameters of the bisector fit to the M 
derived from the sample with temperature profiles. Light grey 
contours describe the slope and normalization of the M — T 
relation for clusters with temperature below 4 keV. In grey we 
show the fit to all the data and in black those excluding the 
groups (M500 < 5 x 10 13 Mq). Contours are drawn at 1,2 and 
3 a confidence level. The innermost contour marks the center. 

Note that it is more straightforward to operate in 
terms of mass in separating the groups, since their tem- 
peratures are strongly affected by even a slight degree of 
preheating (Loewenstein 2000; Tozzi, Scharf & Norman 
2000). 

Restricting ourselves to the systems with temperatures 
above 3 keV, the fit is 



Mi 



r 500 = (4.22I^)10 13 



x hT ew 



The low-temperature end of the M — T relation (sys- 
tems with temperatures below 4.5 keV) yields 



Mi 



500 



(2.45; 



0.39 



)10 1 



1.87 

X K 1 ew 



+0.15 

0.14 
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The above fits and the data are shown in Fig.|| The 
parameters of the fits are compared with the results for 
the HIFLUGCS and simulations of Evrard et al. (1996) 
in Fig.§ 

To study the confidence area of the parameter esti- 
mation for the above fits, we choose a normalization at 
3.5 keV, which makes the determination of the slope less 
dependent on the normalization. We present the values de- 
rived this way in Fig.[|. It can be seen from this figure that 
the total sample is inconsistent with a power law index 
of 1.5 on more than 99.9% level. The high-temperature 
system sample (excluding groups), although revealing a 
flatter index consistent with a value of 1.5, is not strongly 
deviant from the total sample, e.g. the break in the M — T 
relation has only 95% confidence. A steeper slope, derived 
for the low-temperature end of the sample, can still be 
considered as a fluctuation. However, larger errors in the 
parameter determination in the case of inclusion of groups 
are due to the large spread of groups on the M — T rela- 
tion (see Fig.||). So, the meaning of "fluctuation" is that 
a subsample of systems leading to a derivation of the flat 
slope could be drawn from the existing sample at a high 
probability. The origin of the scatter is further discussed 
in Sec.||. 

In many studies the virial radius is suggested as a unit 
of length {e.g. Evrard et al. 1996; Markevitch et al. 1998; 
Cen & Ostriker 1999) to provide a comparison among the 
systems at equal overdensity. For these estimations, the 
luminosity averaged temperature of the cluster is used 
(Markevitch et al. 1998, FDP). Therefore, we provide 
here a relation between rsoo and the luminosity-weighted 
X-ray temperature, derived from the data in Table 0: 



r 50Q = 0.63: 



o.oi 



V kT e , 



where t^qq is in Mpc and kT ew in keV. For a given tem- 
perature this relation implies a 20% smaller value for rsoo 
compared to similar formulae derived from the simulations 
of Evrard et al. (1996). 

3. Polytropic vs Isothermal 

As given in Eq.|l|, the deduced mass depends on the tem- 
perature and the parameters /3 and 7, describing the shape 
of the density and temperature profile. To further eluci- 
date the origin of the behavior of the M — T relation and 
to circumvent the direct dependence of M on T we exam- 
ine in Fig|| the dependence of (3 (correctly (3x 2 / (1 + x 2 ) at 
?"50o) on T. We note four systems with (3 ~ 0.3 that imply 
a steep dependence of beta on the temperature. They are 
mainly responsible for the steep slope of the M — T rela- 
tion in our sample with spatially resolved temperatures. 
In fact, excluding these systems (IC4296, HCG62, HCG51, 
NGC3258), but leaving other groups in, the M—T relation 
is given by 



Excluding two more systems (A262, HCG94), whose 
/3-index is less than 0.5, no further changes in the derived 
parameters of the M — T relation are obtained: 

M 500 = (3.75t°;^)10 13 x kT^ ■ 

All of the three groups in a sample of Nevalainen 
et al. (2000) have a /3-index less or equal to 0.5. The above 
demonstrates the importance of selection effects on the de- 
rived slope of the M — T relation, especially on the group 
scale, originating from apparent scatter in the M — T rela- 
tion on the scales of groups. Observationally, this scatter 
is correlated with flat gas density profiles, which are often 
taken for signs of preheating (Metzler & Evrard 1997). 
In reality, the situation is more complex, since as pointed 
out by Loewenstein (2000), when the effect of the SN feed- 
back is the only one responsible for the observed M — T 
relation, in low-mass systems, a large amount of metals as- 
sociated with SN explosions should be observed and also 
isentropic cores, corresponding to the case of adiabatic col- 
lapse. Although some of the groups do show such features 
(Finoguenov et al. 2001), these are surprisingly not the 
most deviant systems in the M — T relation. 
1 




0.8 1 



8 10 



Fig. 5. Relation between density gradient (defined by 
/3x 2 /(l + x 2 )) and T ew . 

Horner et al. (1999) have found that the dependence 
of f3 on T ew is responsible for the steepening of the de- 
rived M — T in the isothermal assumption, but when the 
temperature profiles are taken into account, the slope be- 
comes 3/2 again. To verify this, we rewrite Eq.jl] in terms 
of overdensity 



M(5,P,j,T(6)) w 2.2 x 10 15 M Q <5 



-l/ 2/3 3/2 7 3/2 T(5) 3/2 )(2) 



^500 — ( 3 -57tg;26 



1 co + 0.05 



where S denotes the chosen overdensity. If the idea of 
Horner et al. (1999) is correct, Eq.|| should contain a 
term counterbalancing the dependence of (3 on T ew , with 
possibilities 7 oc T" w , T(S)/T ew oc T£ w or both. The index 
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a of this (counter-)dependence is —0.26 ± 0.03 according 
to Horner et al. (1999), but weaker dependences of (3 on 
T are also reported {e.g. 0.16 in Vikhlinin et al. 1999). 

As is seen from Figs.|| and 0, the data show two 
trends, which cancel each other: both very hot and very 
cold systems seem to have stronger temperature gradi- 
ents. So, the trend on the high-energy part counterbal- 
ances the weak dependence of (3 on T, while the trend 
in the low-temperature part only reinforces the trend 
observed in (3. Thus, it becomes clear why inclusion of 
groups makes such a drastic difference in the derived 
M — T relation. An overall fit to these figures gives 
7 = (1.080±0.005)r e t - 010±0 - 004) &ndT(6)/T ew = (0.88± 
0.18)Tj-°- 03±0 - 14) . 
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Fig. 6. Relation between 7 and T ew . The data are shown as 
crosses with error bars drawn at the 68% confidence level. 
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In fact, the arguments presented in Horner 
et al. (1999) are not strictly correct. They base their 
conclusion on the fact that 



M ([3) /M (True) 



Fig. 7. Relation between T500 /T ew and T e , 



r 0A £ or an i n di v id ua i cluster, where 
M(j3) is the mass estimate using the /3-model and assum- 
ing isothermality and M(True) is the true mass of the 
cluster. 

However, what is relevant is: 

M((3)/M(True) oc (r/r core ) 0A for the cluster sample, 
fcore oc r(S) oc T 05 where r core is the cluster core 
radius and r(5) is the radius of overdensity 5, 
from which it follows: 

M([3)/M(True) oc (r(5)/r core ) 0A oc (T°- 5 /T - 5 ) - 4 = 
constant. 

This means that if clusters are self similar, the overesti- 
mate of the beta model is the same for all clusters. In view 
of this consideration, we note that there is a very close 
agreement between the parameters of the M — T relation 
determined in our two samples. This agreement, which 
demonstrates that usage of the isothermality assumption 
does not bias the derived parameters of the M — T re- 
lation, can be used to justify the validity of the M — T 
relation, derived for high-redshift samples, where detailed 
temperature measurements are difficult and an assump- 
tion of isothermality is the only choice. 

Following a suggestion of the referee, we examine also 
the effect of assuming isothermality in deriving the total 
mass for the sample with spatially resolved temperature 
measurements. The results are listed in Table 0. The slope 
of the M—T relation, obtained for the total sample, is still 
steeper, while avoiding the systems with flat (3 gives results 
consistent with HIFLUGCS in both slope and normaliza- 
tion. We note that the lowest f3 value for the HIFLUGCS 
sample is 0.44, which again supports the idea of the im- 
portance of the selection of the systems with respect to 
their values of (3. 

Concluding this section, we identify the inclusion of 
systems with low values of (3 as the most important cause 
of the steep slope of the M — T relation, which could be 
overcome by excluding such systems from the sample. As 
we have shown, a flat (3 is not necessarily a unique char- 
acteristic of groups, which most likely implies a different 
importance of preheating in low-mass systems. This is in 
qualitative agreement with the preferential infall scenario, 
proposed by FAD. 

In all the fits, however, the normalization of the 
M — T relation appeares smaller than in the simulations 
of Evrard et al. (1996). In the following, we would like 
to comment on this issue. As pointed out by Nevalaincn 
et al. (2000), the resolution in the simulations of Evrard 
et al. (1996) is insufficient to resolve the cluster cores. One 
can assume, however, that their simulations are correct 
at low overdensities. Since we measure the temperature 
up to the radius of the overdensity chosen for the mass 
calculations, we can directly check this effect by using a 
temperature at the radius of mass determination, (T500), 
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instead of the luminosity-weighted temperatures. Such a 
comparison may also be less affected by preheating, since 
at 6 = 500 no significant variation in the gas fraction is 
seen (Ettori & Fabian 1999, Vikhlinin et al. 1999). Taking 
measured temperatures also avoids many possible effects 
of averaging, discussed in ME00 and is therefore closer to 
the relations predicted for the mass- weighted temperature. 
A fit to the M500 — T500 relation reads as 



M 500 = (3.29±g;i|) 10 13 



1 89+ - 16 
X kT -0.14 



(errors are stated at the 68% confidence level). 

The normalization of M — T500 in the simulations by 
Evrard et al. (1996) is about 9. - 15. x 1O 13 M (consider- 
ing T500 = 0.6 — 0.8Tx), closer to the observed points, but 
still in obvious disagreement (fixing the slope to 1.5 we 
obtain a normalization at 1 keV of TBll^g x 10 13 M Q ). 




500 



Fig. 8. A/500 — T500 relation. Dashed line shows the rescaled 
simulations by Evrard et al. (1996). 



4. Comparison with optical data 

Comparison between the X-ray and the virial mass mea- 
surements, obtained using velocity dispersions, is long 
known to be subject to contradictions, with velocity bias 
considered as the most likely origin. Recently, a new com- 
pilation of optical measurements was made by Girardi 
et al. (1998), where a much more detailed study was car- 
ried out, e.g. different velocity dispersion profiles in clus- 
ters were identified. Since we have measured the masses 
for many clusters in common, we can combine the X-ray 
mass with velocity dispersion measurements to provide a 
comparison with high-resolution dark matter simulations. 

We take high-resolution KCDM (fl M = 1 - A = 0.3, 
as = 1-0; H Q = 70 km s^ 1 Mpc" 1 ) simulations using 
the ART code (Kravtsov, Klypin, Khokhlov 1997), with 
256 3 particles of 1.1 x 10 9 M Q each in a simulation box of 



60/i _1 Mpc (see Gottloeber et al. 1999 for details). The 
particular runs of the ART code we use for comparison are 
characterized by a scale of velocity dispersion with over- 
density of ap oc S om with a residual scatter around the 
best-fit of 20-30%. Within the simulated box, 17 clusters 
have been identified and we build the M — ap relation, 
scaling the measurements done at different overdensities. 
One possible weakness of such a comparison is that in 
the simulations we have studied the dispersion of the dark 
matter, not of the galaxies. 



5, 

V 




400 



600 800 1000 



a p at r g 



km s 1 



Fig. 9. M — ap relation at overdensity of 90. Circles repre- 
sent the results of high-resolution simulations. Black crosses 
represent a correlation of the X-ray mass measurements with 
the velocity dispersion of galaxies at overdensity of 90. Type 
A clusters (Girardi et al. 1998) are shown in black, type B in 
grey and C in light grey (this identification has been made to 
denote different velocity dispersion profiles measured in these 
clusters, which reflect, according to Girardi et al. (1998), that 
type B and C cluster have been formed earlier and are therefore 
more relaxed). 

Optical observations reveal 3 types of clusters, accord- 
ing to their degree of virialization. In comparing the data, 
we have scaled the optical velocity dispersions according 
to the given type of cluster, using identification and scal- 
ing profiles, reported in Girardi et al. (1998). To compare 
our X-ray results with the work of Girardi et al. (1998), 
we rescale our mass estimates to the overdensity of 90, 
using the Navarro, Frenk & White (1996) profile. 

The resulting M — ap relation is presented in Fig.pl 
The outliers on the M — ap relation, A2319, A3558 and 
A400, are already known to have a peculiar structure (Fer- 
etti, Giovannini, Bohringer 1997; Venturi et al. 2000; 
Lloyd-Davies, Ponman, Cannon 2000). The scatter of the 
other points around the simulations is within 30% and 
therefore comparable to the scatter seen in simulations. 
We note, however, that there is a trend of type B and C 
clusters to have a higher velocity dispersion for a given 
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mass, which is consistent with a definition of their type: 
type B and C are older systems, according to modeling 
of Girardi et al. (1998), that should imply a higher for- 
mation redshift. Since such a correction is not introduced 
into the virial mass calculation, this results in a slight bias 
in mass estimates for these clusters, as seen in FigJTT 



a 10 i 



& 10 




X-ray 



500, z 



Fig. 10. Comparison of the X-ray and optical mass estimates. 
Black crosses denote type A clusters in Girardi et al. (1998), 
grey crosses type B and light grey crosses type C. Type A 
clusters show the best agreement with the X-ray mass mea- 
surements. 



5. Implications on the Redshift (Epoch) of cluster 
formation 

The issue of the redshift of cluster formation was strongly 
suggested in studies of Lilje (1992), Kitayama & Suto 
(1996), Voit & Donahue (1998). The effect on the M -T 
relation is such that for a given mass, the systems that 
formed at earlier times should have higher temperatures. 
Since this scenario is in qualitative agreement with trends 
observed by comparing our sample with the simulations of 
Evrard et al. (1996), we decided to estimate whether the 
shift observed in the M — T relation could be explained 
just by this scenario. 

To do this, we invert the problem, i.e. use the X-ray 
mass and temperature measurements and the M — T rela- 
tion obtained in the simulations to derive the distribution 
of redshifts of cluster formation for our sample. In more 
detail, for a measured mass of the system, we compare 
the measured temperature with the value obtained from 
the simulation for the measured mass and attribute the 
difference to the redshift of cluster formation, which for 
fi M = 1 is simply T observed = T simulated x (1 + z f ). To 
be able to do this, one should make sure that the simu- 
lated relation explicitly assumes that the clusters form at 
the redshift of observation. Fortunately, the simulations of 



Evrard et al. (1996) have this assumption, which is quite 
logical for J7m = 1, used for most of their runs (Metzler, 
private communication). 

In Fig.|ll| we illustrate the effect of the redshift of clus- 
ter formation, by comparison with the model of Lacey & 
Cole (1993) following the formulae presented in Balogh 
et al. (1999) for flu = 0.3. We subdivide our sample 
into 3 parts with masses in the 0.1 — 0.8, 0.8 — 3 and 
3 — 15 x 10 14 Mq range. From the Figure one can see that 
the theoretical prediction varies significantly between the 
three subsets and generally matches the trends seen in the 
corresponding subset. While a more detailed comparison 
should await an M — T relation simulated with much bet- 
ter resolution, we can state already that a steeper slope of 
the observed M — T relation implies that lower-mass sys- 
tems form preferentially at earlier times, compared to rich 
clusters. The formation redshift distribution for groups is 
wider than the model prediction, which can be taken as 
a sign of the importance of SN preheating, but a more 
detailed study is needed to verify this suggestion. 



1 

* 
N 

2 



T3 



I 1 1 1 1 I 

M ~ 5x10" M„ 



1TI Mil l 



M ~ 1x10'" M„ 




i i i |~ t~ i ~ s t'H- 



M ~ 5x10" «„ 



0.5 



1.5 



Fig. 11. Redshifts of cluster formation, deduced from a re- 
quirement for measurements to match the simulations of 
Evrard et al. (1996). The black solid histogram shows the cal- 
culation for the sample with temperature profiles with dashed 
histograms indicating the number count errors at the 68% 
confidence level. We present the comparison with the model 
of Lacey & Cole (1993) following the formulae presented in 
Balogh et al. (1999) for Qm = 0.3 and three mass ranges in 
our sample with typical masses of 5 x 10 13 Mq, 1 x 10 14 Mq 
and 5 x 10 14 Mq (grey lines). 

We can already identify one source of bias in the for- 
mation redshift distribution of the clusters in the present 
study. Systems that collapsed recently bear indications of 
recent merger activity and therefore are selectively re- 
moved from the sample aimed to determine the mass 
using the assumption of hydrostatic equilibrium. There- 
fore, the first bin in the derived Zf distribution should be 
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disregarded in FigJT^. Flux-limited samples, such as HI- 
FLUGCS, arc in principle also biased in this regard, since 
older systems are expected to be brighter for a given mass. 
The lack of clusters lying close to the simulations on the 
M — T plane is also caused by influence of a central region 
on the determination of the luminosity averaged temper- 
ature, since formation of the central part of a cluster is 
slightly shifted towards higher formation epochs. 



6. Conclusions 

We have studied the M — T relation using a sample of 
clusters with resolved temperature profiles and also using 
HIFLUGCS. Our findings are 



Using spatially resolved temperatures in mass esti- 
mates results in similar parameters of the M — T re- 
lation to the analysis assuming isothermality. The ob- 
served slope of the M — T relation is steeper than the 
predictions of the self-similar relations and the mea- 
sured normalization is two times lower than that ob- 
tained in the simulations of Evrard et al. (1996). 
Significant scatter of the points on the M — T relation 
was found for groups. This is a likely source of the 
disagreement found between the studies of Nevalainen 
et al. (2000) and Horner et al. (1999). 
We find that the derived slope of the M — T relation 
depends strongly on the (3 dependence on temperature 
and is not counterbalanced by inclusion of temperature 
gradients, as proposed by Horner et al. (1999). This 
provides an explanation for similar parameters of the 
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Table 2. Mass determinations using spatially resolved temperatures. 



Name 


z 


kT 


fcT 50 o 


M 500 


7"500 


a 


f' 'core 


7 


Rout 






keV 


keV 


10 14 M e 


Mpc 




Mpc 




Mpc 


A2029 


0.077 


9.10 ±0.66 


7.87 ±0.86 


12.71 ±2.44 


2.06 ±0.27 


0.68 


0.28 


1.11 ±0.08 


2.82 


A401 


0.074 


8.00 ±0.24 


5.34 ± 0.47 


9.40 ± 1.32 


1.70 ±0.16 


0.63 


0.27 


1.20 ±0.07 


1.82 


A3266 


0.055 


8.00 ±0.30 


5.61 ±0.98 


10.36 ±2.81 


1.93 ±0.36 


0.74 


0.50 


1.29 ± 0.14 


1.57 


A 1795 


0.062 


7.80 ± 0.60 


5.77 ± 1.39 


11.66 ±4.72 


2.00 ±0.56 


0.83 


0.39 


1.19 ±0.19 


1.56 


A2256 


0.058 


7.30 ±0.30 


4.64 ± 0.62 


8.38 ± 1.80 


1.79 ±0.27 


0.82 


0.52 


1.26 ±0.11 


1.83 


A3571 


0.040 


6.90 ±0.12 


5.12 ± 1.93 


6.91 ±2.75 


1.68 ±0.46 


0.69 


0.27 


1.12 ±0.08 


2.28 


A1651 


0.085 


6.10 ±0.24 


5.05 ± 1.09 


6.77 ± 1.85 


1.67 ±0.32 


0.70 


0.26 


1.11 ±0.11 


1.62 


A119 


0.044 


5.60 ±0.18 


4.96 ±0.39 


5.34 ±0.95 


1.54 ±0.19 


0.66 


0.48 


1.09 ± 0.11 


1.28 


A3558 


0.048 


5.50 ±0.24 


4.44 ± 0.93 


4.05 ± 1.01 


1.41 ±0.24 


0.55 


0.19 


1.14 ±0.10 


1.39 


A2199 


0.030 


4.80 ±0.12 


4.19 ±0.29 


3.95 ±0.41 


1.40 ±0.10 


0.64 


0.14 


1.00 ± 0.05 


1.25 


A496 


0.033 


4.70 ±0.12 


3.84 ± 0.46 


4.70 ±0.67 


1.48 ±0.15 


0.70 


0.25 


1.15 ±0.05 


1.32 


A4059 


0.048 


4.40 ±0.18 


3.94 ± 0.77 


3.87 ± 1.20 


1.39 ±0.30 


0.67 


0.22 


1.02 ±0.17 


0.92 


A3112 


0.075 


4.34 ±0.30 


3.74 ±0.41 


3.54 ± 0.45 


1.35 ±0.12 


0.63 


0.12 


1.06 ± 0.04 


1.74 


Hydra 


0.057 


4.30 ±0.24 


3.85 ±0.51 


3.89 ±0.71 


1.39 ±0.18 


0.66 


0.12 


1.05 ± 0.09 


1.35 


A2063 


0.035 


3.86 ±0.14 


3.55 ±0.46 


3.68 ± 0.54 


1.36 ±0.14 


0.69 


0.22 


1.07 ±0.04 


0.93 


MKW3S 


0.045 


3.79 ±0.32 


2.18 ±0.26 


2.57 ±0.43 


1.21 ±0.14 


0.71 


0.30 


1.37 ±0.12 


1.16 


A2657 


0.040 


3.70 ±0.18 


2.30 ±0.69 


2.93 ± 1.02 


1.27 ±0.31 


0.76 


0.37 


1.36 ±0.19 


1.17 


AWM7 


0.017 


3.50 ±0.12 


3.54 ±0.27 


2.33 ±0.21 


1.17 ±0.07 


0.53 


0.10 


1.01 ±0.02 


0.75 


A2052 


0.035 


3.46 ± 0.24 


3.54 ±0.79 


3.02 ± 0.96 


1.28 ±0.28 


0.64 


0.10 


0.99 ±0.17 


0.46 


HCG94 


0.042 


3.45 ±0.31 


2.68 ±0.36 


1.67 ±0.32 


1.05 ±0.14 


0.48 


0.08 


1.16 ± 0.14 


1.08 


2A0335 


0.035 


3.40 ± 0.24 


3.95 ± 1.20 


3.44 ± 1.13 


1.34 ±0.30 


0.65 


0.08 


0.95 ± 0.08 


0.92 
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0.028 


3.31 ±0.16 


3.55 ±0.28 


2.55 ±0.50 


1.21 ±0.16 


0.61 


0.16 


0.94 ±0.14 


0.38 
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3.14 ±0.06 


1.99 ±0.80 


1.98 ±0.86 


1.11 ±0.33 


0.70 


0.16 


1.23 ±0.17 


0.31 
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2.60 ±0.81 


1.22 ±0.26 


0.69 
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1.35 ±0.27 


0.83 
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2.92 ±0.26 


3.10 ±0.57 


1.77 ±0.94 


1.07 ±0.39 


0.52 
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0.97 ± 0.45 


0.71 
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2.63 ±0.48 


1.67 ±0.43 


1.05 ±0.19 
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2.66 ±0.27 


2.28 ± 0.44 


1.16 ±0.16 


0.69 
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1.05 ± 0.14 


0.72 
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0.028 


2.29 ±0.16 


1.91 ±0.18 


1.03 ±0.17 


0.89 ±0.10 


0.51 


0.12 


1.14 ±0.16 


0.75 


A262 
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2.26 ±0.14 


2.13 ±0.53 


0.93 ±0.27 


0.86 ±0.17 
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1.04 ±0.17 
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1.29 ±0.42 
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0.57 
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1.55 ±0.43 


0.91 ±0.26 


0.86 ±0.17 
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0.18 


1.06 ± 0.02 


0.66 


N6329 


0.028 


1.45 ±0.08 


1.44 ± 0.40 


0.68 ±0.19 


0.78 ±0.15 


0.53 


0.12 


1.10 ±0.03 


0.64 


HCG51 


0.026 


1.34 ±0.07 


1.44 ± 0.43 


0.22 ±0.07 


0.54 ±0.11 


0.30 


0.08 


0.93 ±0.10 


0.69 


N5044 


0.009 


1.23 ±0.04 


0.87 ±0.07 


0.36 ±0.07 


0.63 ±0.08 


0.51 


0.01 


1.22 ±0.30 


0.28 


HCG62 


0.014 


1.11 ±0.05 


0.84 ±0.20 


0.19 ±0.05 


0.51 ±0.09 


0.31 


0.02 


1.36 ±0.19 


0.58 


N4325 


0.026 


0.98 ± 0.04 


0.86 ±0.45 


0.40 ± 0.25 


0.65 ±0.28 


0.59 


0.01 


1.14 ±0.56 


0.60 


IC4296 


0.013 


0.95 ±0.09 


0.73 ±0.35 


0.15 ±0.07 


0.47 ±0.16 


0.31 


0.06 


1.35 ±0.15 


0.32 


N5129 


0.023 


0.75 ±0.05 


0.45 ±0.09 


0.23 ±0.05 


0.54 ±0.08 


0.60 


0.10 


1.51 ±0.04 


0.58 



